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1 .  Introduction 

Axisymmetric  rotating  flows  ov?r  a  rotating  disk  of  infinite  di¬ 
mensions  have  been  studied  in  [ 1  through  8].  Recently  the  special 
problems  of  von  Kirm^n  and  Bodewadt  have  been  reinvestigated  in  [9 
and  12],  Adjusted  local  boundary  layer  approximations  of  first  order, 
which  were  introduced  in  [10  and  11],  have  been  obtained  for  both 
problems.  The  heat  transfer  from  a  rotating  disk  into  a  fluid 
asymptotically  at  rest,  which  is  computed  in  [9],  agrees  well  with 
empirical  results. 

In  the  present  paper  the  local  boundary  layer  approximations  of 
first  order  derived  In  [12]  are  generalized  to  axisymmetric  motions 
which  rotate  over  a  rotating  disk  of  infinite  dimensions.  Numerical 
results  are  computed  and  discussed  for  a  variety  of  Reynolds  numbers 
and  for  cases  in  which  the  disk  is  rotating  in  the  same  sense  and  In 
the  opposite  sense  as  the  fluid  far  away  from  the  disk.  In  particular, 
critical  Reynolds  numbers  for  steady  laminar  motions  which  are  attached 
to  the  surface  of  the  disk  are  computed  and  displayed. 

2 .  Construction  of  Rotating  glows 

A  viscous  incompressible  fluid  may  be  contained  in  a  circular 
cylinder  which  is  normal  to  a  disk  of  equal  diameter  (see  Fig.  1). 

If  the  cylinder  is  rotating  at  a  constant  angular  velocity  0), 

then  the  fluid  far  away  from  the  disk  will  essentially  assume  the 
same  angular  velocity  u>».  Near  the  disk  the  fluid  rotates  with  the 
angular  velocity  0)  of  the  disk.  Furthermore,  due  to  the  friction 

forces  at  the  disk, fluid  will  be  ejected  or  drawn  through  a  slit  of 
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negligible  size  which  may  be  provided  between  the  cylinder  and  the  disk. 
Hence,  the  resulting  rotation  of  the  fluid  becomes  superposed  by  a 
secondary  flow  as  sketched  in  Figs.  la.  through  1  c.. 

In  order  to  investigate  exclusively  the  effects  of  the  friction 
forces  at  the  disk  on  the  rotating  fluid,  the  surface  of  the  cylinder 
wall  may  be  assumed  to  glide  freely  with  the  generated  axial  flow.  For 
the  sake  of  simplicity,  only  the  flow  model  of  infinite  diameter  and 
height  is  considered  which  is  dynamically  similar  to  one  of  the  finite 
models  described.  As  was  pointed  out  in  [12]  the  infinite  flow  model 
is  specified  by  the  constant  kinematic  viscosity  v  of  the  fluid,  by 
the  two  angular  velocities  cc0  and  coo,,  and  by  the  characteristic  length 
X(>  0)  of  the  motions  considered.  The  latter  parameter  (see  [12])  may 
be  defined  as  the  geometric  mean  of  the  axial  thickness  a  and  twice  the 
axial  radius  of  curvature  1/b  of  a  properly  defined  boundary  layer  along 
the  solid  disk  (see  Eqs.  (10)  and  (12)).  Accordingly,  a  finite  flow 
model  and  an  infinite  model  are  considered  geometrically  similar  if 
they  agree  in  their  characteristic  lengths. 

For  the  mathematical  description  of  the  flow  a  coaxial  cylindrical 
coordinate  system  (r,cp,z)  may  be  introduced  so  that  z  =  0  designates  the 
solid  surface  of  the  disk.  If  (u,v,w)  denotes  the  corresponding  velocity 
vector,  p  the  variable  pressure,  and  p  the  constant  density,  then  the 
motion  is  governed  by  the  Navier-Stokes  equations: 

uur  +  wuz  -  ^  |  Pr  +  v[urr  +  (?)r  +  uzz]  (D 
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uvr  +  wvz  +  if 


uwr  4-  wwz 


+  Vr„._  +  M  +  vi 

!  \l/i  J-“J 

-  -  i  pz  +  v[wrr  +  £  wr  +  w2Z] 


(3) 


(ru)r  +  (rw)z  -  0  .  (4) 

In  order  to  introduce  dimensionless  quantities  it  is  useful  to 
define 

-  Wq  -  lUoo  and  *  ■  —  (5) 

respectively  as  the  "characteristic  angular  velocity"  and  the  "charac¬ 
teristic  velocity  ratio"  of  the  motion  considered.  In  contrast  to  other 
definitions  (see  [1,  6,  8])  the  expressions  (5)  appear  to  be  advantageous 
because  they  become  meaningless  (Q  “  0,  k  *=  °°)  only  for  the  trivial  ro¬ 
tation  which  is  defined  by  uu0  *  x®. 

In  the  dimensionless  depend  it  variables 


u  ■  QrU,  v  ■  QrV,  w  -  v<JW,  — 


P  _  fl2  r2  ..a 


h  TP 


(6) 


the  solution  of  the  Navier-Stokes  equations  under  consideration  is 
determined  by  the  boundary  data  (see  [12]) 

z  -  0,  |r|  <  »:  U  -  0,  V  -  1  +  h,  W  -  0  ,  (7) 

z  -  ®,  ]  r  |  s:  ®  I 

f :U  -  0,  v  -  k  ,  W  -  1,  (8) 

z  >  0,  | r ]  *  ®  * 

where  w^,  is  a  certain  constant  which  must  be  determined  simultaneously 
with  the  integral. 

Following  derivations  presented  in  [12]  it  is  helpful  to  Introduce 
the  limiting  line  of  the  boundary  layer  along  the  disk 
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z  -  6(r)  -  a  -  br2  -I-  cr4  +  ...  (r  <  ru)  (9) 

by  the  e-condition 

z  *  6(r)  :  jv  -  h]  ■  e  (e  >  0)  .  (10) 

A  local  boundary  layer  approximation  of  first  order,  which  is  adjusted 
to  the  curvature  of  the  boundary  layer  around  the  axis  of  rotation  r  =  0, 
is  of  the  form 

*  2  a  ^ 

U  =  -  GOO,  V  =  V  (T[)  ,  w  =  —  1G(7]),  P  =  -Z-  ^2H(T),  (11) 

provided  the  motion  is  steady,  laminar,  and  attached  to  the  disk.  In 
these  equations  the  dot  denotes  the  derivative  with  respect  to  the  single 
dimensionless  variable 

T1  “  M  77-  =  - -  ,  £  =  \/a/b  (12) 

6(r)  1  -  f  +  ... 

where  £  is  by  definition  the  characteristic  length  of  the  motion.  After 
f  12]  the  first  order  approximation  (11)  is  the  integral  of  the  differential 
equations 

G  +  8T|G  =*  -  R(;t2  +  G2  -  2GG  -  V2  -  4T[H)  (13) 

V  +  871V  -  -  2R(GV  -  GV)  (14) 

H  -  2GG  -  i  (G  +  4T|G)  (15) 

K 

which  satisfies  the  boundary  conditions 

T)  =  0:  G  0,  G  -  0,  V  =  1  +  k  (16) 

T|  "  cn:  G  -  Gm,  V  =  h ,  H  =  G2  .  (17) 

The  dimensionless  parameter 

R  =  7  £2  (18) 

can  be  interpreted  as  the  Reynolds  number  of  the  motion.  It  may  be 
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mentioned  that  by  this  definition  the  Reynolds  number  assumes  posicive 
and  negative  values  according  as  fi  is  positive  or  negative  (see  Eq.  (5)). 
This  conventic  appears  to  be  useful  since  the  axial  flow  is  expected  to 
be  of  stagnation  type  (directed  against  the  disk)  for  fi  =  uj0  -  ujo,  >  0 
(see  [13]).  For  0  -  uu0  -  oj®  <  0  the  axial  flow  will  generally  be  of 
wake  type  (directed  away  from  the  disk)  with  the  few  exceptions  dis¬ 
played  in  Figs.  1  a.  through  Id..  It  may  also  be  noted  that  by  con¬ 
vention  (5)  the  characteristic  velocity  rvtio  h  possesses  the  so.  e 
sign  as  fi  and  R.  In  particular,  the  flows  of  von  Kirm^n  and  Bodewadt, 
which  have  been  investigated  in  [9  and  12],  are  characterized  by  the 
parameters  (R  >  0,  x  «  0)  and  (R  <  0,  h  -  -  1),  respectively. 

3 .  Stability  of  Rotating  Flows 

Axisymmetric  rotating  flows  over  a  rotating  disk  have  been  computed 
for  various  Reynolds  numbers  R  and  velocity  ratios  x  (Tables  1-12).  Se¬ 
lected  examples  are  displayed  in  Figs.  2  through  14,  As  in  von  K^rmin's 
case  h  ■  0  (see  [12])  numerical  calculations  of  rotating  flows  with 
h  *  0  (i.e.,  uj0  >  (jjot  >  0)  indicated  no  symptoms  of  instability  for  any 
Reynolds  number  R(>  0).  The  secondary  flow  is  of  pure  stagnation  type 
(see  Fig.  la)  and  all  velocity  profiles  are  free  of  inflection  points 
which  characterize  unstable  motions. 

For  velocity  ratios  k  <  0  (i.e.,  uj0  <  (A»(>  0))  all  rotating  flows 
become  unstable  beyond  a  certain  "critical  Reynolds  number"  Rc(h)(<  0) 
which  depends  on  the  value  of  h.  Fig.  2  shows  the  corresponding 
critical  curve  which  has  been  approximately  determined  at  several  points. 
These  numerical  calculations  demonstrated  distinctly  two  different  kinds 
of  instabilities  below  and  above  x  «  -  1/2  (i.e.,  (Da,  «  -  uq,  ) . 
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For  x  <  -  1/2  (i.e.,  ud^  >  |u>0|)  the  instability  is  of  the  same  sort 
as  in  Bodewadt's  case  h  =  -  1  (see  [12]).  The  secondary  flow  is  of  pure 
wake  type  (see  Fig.  lc.)  which  begins  to  oscillate  in  the  upper  portion 
of  the  boundary  layer  when  the  Reynolds  number  R  exceeds  its  critical 
value  Rc.  This  experience  indicates  that  the  oscillating  solutions, 
which  were  found  in  [1,  6,  8]  for  x  £  -  1  (i.e.,  u)ro  >  u)0  0) ,  should 

also  exist  for  the  x-range  -  1/2  £  x  <  -  1.  Unfortunately,  Rogers  and 
Lance  [6]  inspected  the  entire  x-range  with  the  exception  of  this  inter¬ 
val.  However,  as  was  pointed  out  in  [12  and  13]  the  nonexistence  of  a 
proper  solution  to  the  boundary  value  problem  (13)  through  (17)  is  an 
indication  that  in  reality  the  motion  is  separated  from  the  surface  of 
the  disk  (see  the  tea-cup  experiment  described  in  [12])  and  violates 
the  hypothesis  (11). 

For  the  x-range  -  l/2<x<0  (i.e.,  0  <  Wc  <  -  )  the  instability 

beyond  a  critical  Reynolds  number  Rc(x)  arises  in  the  lower  portion  of 
the  boundary  layer  through  the  occurrence  of  two  inflection  points  in 
the  tangential  velocity  profile  (see  [13]  ).  The  solutions  of  the  dif¬ 
ferential  equations  (13),  (14),  and  (15)  appear  to  be  unbounded  at  infinity 
and  it  is  impossible  to  satisfy  the  boundary  conditions  (17).  In  the 
x-range  considered  the  secondary  flows  are  either  of  the  "cell  type" 
sketched  in  Figs.  I'd,  and  lc.  or  of,  the  "distorted  stagnation  type" 
plotted  in  Fig.  lb.,  provided  the  absolute  value  of  the  Reynolds 
number  R  is  sufficiently  large.  When  the  Reynolds  number  R  exceeds 
the  critical  value  Rc(x),  the  whole  flow  field  seems  to  be  affected 
by  the  friction  forces  at  the  disk.  This  explains  the  nonexistence 
of  a  proper  solution  to  the  problems  considered  in  [1,  6,  8  ]  for 
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the  same  n-range  0  <  k  £  -  1/2.  Indeed,  these  problems  are  specified 
among  the  present  problems  by  infinite  Reynolds  numbers.  Since  k  -  -  1/6 
(i.e.,  u)0  •  -  5u)oo(<  0))  marks  roughly  the  "limiting  cell  flow"  sketched 
in  Fig.  lc.  for  large  Reynolds  numbers,  the  instability  will  be  most 
severe  for  the  cell-type  range  -1/2shs-1/6  and  less  severe  for 
the  distorted  stagnation- type  range  -1/6<h<0.  In  fact,  this  con¬ 
firms  and  explains  the  numerical  experience  of  Rogers  and  Lance  as  de¬ 
scribed  in  [6], 

For  practical  applications  it  is  important  to  observe,  that  accord¬ 
ing  to  the  curve  of  critical  Reynolds  numbers  presented  in  Fig.  2,  the 
instability  of  rotating  motions  increases  with  increasing  velocity 
ratios  (-  k)  >  0,  i.e.,  with  decreasing  characteristic  angular  velocity 
(-  fi)  ■  ala,  -  (u0  >  0,  This  significant  phenomenon  is  physically  plausible. 
Indeed,  if  a  fluid  rotates  at  an  angular  velocity  u)ro(>  0)  over  a  disk, 
the  tendency  of  flow  separation  will  be  supported  by  a  slow  (0  <  <  tuo,) 

rotation  of  the  disk  in  the  same  direction.  The  rotation  of  the  disk 
in  the  direction  of  the  rotating  fluid  acts  similarly  to  an  injection 
of  fluid  through  the  disk  which  decreases  the  pressure  variation  along 
the  axis  of  symmetry  (see  Fig.  6a.)  and,  thus,  increases  the  danger  of 
flow  separation.  On  the  other  side,  the  tendency  of  flow  separation  can 
be  greatly  reduced  by  rotating  the  disk  in  opposite  direction  of  the 
rotating  fluid.  In  this  case  the  rotation  of  the  disk  acts  like  a 
suction  of  fluid  through  the  disk  which  increases  the  pressure  vari¬ 
ation  along  the  axis  of  symmetry  (see  Fig.  6a.)  and,  hence,  decreases 
the  danger  of  flow  separation. 
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The  phenomenon  of  flow  separation  can  eventually  be  eliminated  by 

rotating  the  disk  sufficiently  fast  in  opposite  sense  of  the  rotation 
of  the  fluid.  However,  such  a  rotation  produces  cell  flows  or  even 
distorted  stagnation  flows  (see  Figs.  Id.,  1c.,  and  lb.)  which  reach  a 
different  kind  of  critical  point.  When  the  Reynolds  number  approaches 
its  critical  value,  the  friction  forces  at  the  surface  affect  the  en¬ 
tire  flow  field  and  prevent  a  proper  solution  from  existing. 

Flow  Patterns  of  Rotating  Motions 

In  Figs.  3a.  through  6b.  the  velocity  profiles  V,  G  =  w^W/Zf 

• 

G  ■  -  U  and  the  pressure  P  have  been  plotted  for  various  velocity 
ratios  x,  which  have  been  computed  at  the  constant  Reynolds  number 
| R |  ■  1.  While  the  tangential  velocity  profiles  V  are  essentially 
congruent  for  all  velocity  ratios  x,  the  secondary  velocity  profiles 
G  and  G  display  significant  changes  especially  in  the  most  interesting 
x-range  -  1/3  £  x  s  0.  From  the  secondary  velocity  profiles  G  and  G 
one  arrives  easily  at  the  following  five  essentially  different  second¬ 
ary  flow  patterns. 

a.  Pure  Stagnation  Flows:  At  I R j  =1  the  secondary  flow  is  of 
pure  stagnation  type  within  the  x-range  -  .035  £  x  <  +  ®  (see  Fig.  la.). 
The  pressure  displays  the  same  nonmonotonic  behavior  as  in  von  Kirmin's 
case  x  ■  0  (see  [12]  and  Figs.  6a.  and  6b.)  which  represents  the  proto¬ 
type  of  these  motions.  It  may  be  mentioned  that  the  lower  limit 

x  -  -  .035,  which  has  been  approximately  determined, increases  to 
x  ■  0  (u^,  >  ujco  >  0)  when  the  Reynolds  number  ! R I  increases  to  infinity. 

b.  Distorted  Stagnation  Flows;  For  the  x-range  -  .2056  <  x  <  -  .035 
and  R  -  -  1  the  secondary  flows  are  of  the  distorted  stagnation  type 
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sketched  in  Fig.  lb.  In  this  x-range  the  pressure  is  subjected  to  its 

greatest  change  (see  Fig.  6b.).  When  k  varies  from  the  upper  limit 

(h  ■  -  .035)  to  the  lower  limit  (x  ■  -  .2056)  the  pressure  changes 

rapidly  from  the  stagnation  type  to  the  wake  type  and  maintains  that 

feature  for  all  lower  values  of  x.  In  Figs.  7.  through  10.  a  cross- 

section  of  distorted  stagnation  flows  is  shown  for  various  Reynolds 

numbers  R  at  the  fixed  velocity  ratio  x  =  -  .1.  One  can  see  that  for 

the  velocity  ratio  x  *  -  .1  the  stagnation  character  of  the  secondary 

flow  remains  unchanged  for  all  noncritical  Reynolds  numbers.  In  fact, 

the  computations  indicated  that  the  x-range  of  distorted  stagnation 

flows  changes  to  approximately  the  interval  -  l/6<x<0  (0<  Su^  <  -  ) 

if  the  Reynolds  numbers  R  approach  their  critical  values  Rc(x). 

c.  Limiting  Cell  Flows:  For  R  ■  -  1  and  x  a;  -  .2056  the  secondary 
motion,  which  is  always  generated  by  the  rotating  fluid  at  large  dis¬ 
tances  from  the  disk,  vanishes  completely  and  a  simple  cell  flow 
develops  as  sketched  in  Fig.  lc.  The  pressure  assumes  essentially 

the  nonmonotonic  features  of  wake  flows  which  were  pointed  out  for 
Bodewadt  flows  (x  ■  -  1)  in  [12],  When  the  Reynolds  number  |R| 
increases,  the  velocity  ratio  x  of  limiting  cell  flows  approaches 
the  value  x  cs  -  1/6  (5ajo  ss  -  uu0 ) . 

d.  Cell  Flows;  Within  the  x-range  -  .3271  <  x  <  -  .2056  the 
friction  forces  at  the  rotating  disk  are  still  strong  enough  to  produce 
a  kind  of  degenerated  stagnation  flow  near  the  disk.  However,  far 
away  from  the  disk  the  rotation  of  the  fluid  dominates  over  the  effects 
of  the  rotation  of  the  disk  and  produces  a  wake  flow  so  that  a  cell  flow 
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develops  as  sketched  in  Fig.  Id.  Hence,  while  the  pressure  undergoes 

its  greatest  change  in  the  x-range  of  distorted  stagnation  flows,  the 

secondary  motion  suffers  its  greatest  change  in  the  x-range  of  cell 

flows.  The  graphs  in  Figs.  11.  through  14.  display  a  cross-section  of 

cell  flows  for  various  Reynolds  numbers  R  at  the  fixed  velocity  ratio 

k  *  -  .22.  As  these  graphs  show  the  cell  type  remains  unchanged  for 

all  noncritical  Reynolds  numbers  provided  x  =  -  .22.  Computations 

indicated  that  the  x-range  of  cell  flows  increases  to  approximately 

the  interval  -  1 /2  <  k  <  -  1/6  (0  <-ub/5  <  uuro  <  -  uu  )  when  the 

Reynolds  number  increases  in  its  absolute  value  'Rl. 

e.  Pure  Wake  Flows:  In  the  x-range  -  ®  <  h  <  -  .3271  at  R  =  -  1 

the  rotation  of  the  fluid  far  away  from  the  disk  determines  the  entire 

secondary  flow.  A  pure  wake  flow  develops  as  sketched  in  Fig.  lc. 

The  prototype  of  these  secondary  flows  is  the  wake  flow  which  was 

described  in  [12]  in  Bodewadt's  case  x  =  -  1 .  When  the  Reynolds 

number  1 R I  increases  toward  its  critical  value  Rc(h),  the  x-range 

of  pure  wake  flows  seems  to  approach  the  semi-infinite  interval 

—  ™<H<—  1/2  (uOoo  >  J  U^-j  1 )  • 
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TABLE  4:  PRESSURE  DISTRIBUTION  P  (|R|  -  1) 


TABLE  5 


TANGENTIAL  VELOCITY  V  FOR  k  -  .1 


71 

R  -  -  1 

R  -  -  10 

R  =  -  50 

R  =  -  100 

.00 

.9000 

.9000 

.9000 

.9000 

.9000 

.05 

.7875 

.7874 

.7780 

.7073 

.6435 

.10 

.6773 

.6771 

.6596 

.5352 

.4336 

.15 

.5714 

.5711 

.5476 

.3906 

.2760 

.20 

.4716 

.4713 

.4442 

.2732 

.1613 

.25 

.3795 

.3792 

.3504 

.1796 

.7882-1 

.30 

.2961 

.2958 

.2672 

.1062 

.2020-1 

.35 

.2222 

.2219 

.1946 

.4956-1 

-.2089-1 

.40 

.1579 

.1576 

.1326 

.6540-2 

-.4912-1 

.45 

.1031 

.1028 

.8070-1 

-.2552-1 

-.6807-1 

.50 

.5730-1 

.5706-1 

.3802-1 

-.4892-1 

-.8046-1 

.55 

.1980-1 

.1959-1 

.3635-2 

-.6565-1 

-.8834-1 

.60 

-.1031-1 

,1048-1 

-.2352-1 

-.7735-1 

-.9322-1 

.65 

-.3401-1 

-.3414-1 

-.4454-1 

-.8536-1 

-.9616-1 

.70 

-.5228-1 

-.5239-1 

-c 604 9-1 

-.9072-1 

-.9788-1 

.75 

-.6610-1 

-.6618-1 

-.7235-1 

-.9423-1 

-.9886-1 

.80 

-.7635-1 

-.7641-1 

-.8100-1 

-.9648-1 

-.9940-1 

.85 

-.8379-1 

-.8383-1 

-.8718-1 

-.9790-1 

-.9969-1 

.90 

-.8909-1 

-.8912-1 

-.9150-1 

-.9877-1 

-.9985-1 

.95 

-.9279-1 

-.9281-1 

-.9447-1 

-.9929-1 

-.9993-1 

1.00 

-.9532-1 

-.9534-1 

-.9647-1 

-.9960-1 

-.9997-1 

1.10 

-.9814-1 

-.9814-1 

-.9864-1 

-.9988-1 

-.9999-1 

1.20 

-.9931-1 

-.9931-1 

-.9951-1 

-.9997-1 

-.9999-1 

1.30 

-.9976-1 

-.9976-1 

-.9984-1 

-.9999-1 

1.40 

-.9992-1 

-.9993-1 

-.9995-1 

-.9999-1 

TABLE  6 


AXIAL  VELOCITY  G  -  VcoW/20e  FOR  h  -  -  .1 


T) 

R  -  -  .1 

R  -  -  1 

R  =  -  10 

R  =  -  50 

R  -  -  100 

.00 

0 

0 

0 

0 

0 

.05 

.9346-5 

.9334-4 

.8624-3 

.2577-2 

.3519-2 

.10 

.3171-4 

.3167-3 

.2893-2 

.7914-2 

.9911-2 

.15 

.6009-4 

.6001-3 

.5424-2 

.1361-1 

.1570-1 

.20 

.8947-4 

.8935-3 

.7992-2 

.1848-1 

.1980-1 

.25 

.1166-3 

.1165-2 

.1031-1 

.2213-1 

.2223-1 

.30 

.1398-3 

.1396-2 

.1225-1 

.2457-1 

.2342-1 

.35 

.1583-3 

.1580-2 

.1376-1 

.2605-1 

.2384-1 

.40 

.1722-3 

.1719-2 

.1486-1 

.2683-1 

.2384-1 

.45 

.1821-3 

.1818-2 

.1562-1 

.2715-1 

.2366-1 

.50 

.1887-3 

.1884-2 

.1611-1 

.2722-1 

.2344-1 

.55 

.1929-3 

.1925-2 

.1640-1 

.2715-1 

.2324-1 

.60 

.1953-3 

.1950-2 

.1656-1 

.2704-1 

.2309-1 

.65 

.1966-3 

.1962-2 

.1663-1 

.2692-1 

.2297-1 

.70 

.1971-3 

.1967-2 

.1665-1 

.2681-1 

.2290-1 

.75 

.1972-3 

.1968-2 

.1664-1 

.2673-1 

.2285-1 

.80 

.1972-3 

.1968-2 

.1663-1 

.2668-1 

.2282-1 

.85 

.1970-3 

.1966-2 

.1661-1 

.2664-1 

.2280-1 

.90 

.1968-3 

.1964-2 

.1660-1 

.*661 -1 

.2279-1 

.95 

.1966-3 

.1963-2 

.1658-1 

.2659-1 

.2278-1 

1.00 

.1965-3 

.1962-2 

.1657-1 

.2659-1 

.2278-1 

1.10 

.1964-3 

.1960-2 

.1656-1 

.2658-1 

.2278-1 

1.20 

.1963-3 

.1959-2 

.1656-1 

.2657-1 

.2277-1 

1.30 

.1963-3 

.1959-2 

.1656-1 

.2657-1 

1.40 

.1963-3 

.1959-2 

.1655-1 

.2657-1 

TABLE  7 


RADIAL  VELOCITY  G  «  -  U  FOR  x  -  -  .1 


11  R--.1  R  -  -  1  R  =  -  10  R  =  -  50  R  =  -  100 


.00 

0 

0 

0 

0 

0 

.05 

.3438-3 

.3434-2 

.3153-1 

.8963-1 

.1162 

.10 

.5273-3 

.5267-2 

.4748-1 

.1160 

.1285 

.15 

.5912-3 

.5904-2 

.5221-1 

.1081 

.1000 

.20 

.5735-3 

.5725-2 

.4959-1 

.8558-1 

.6417-1 

.25 

.5068-3 

.5058-2 

.4285-1 

.6036-1 

.3464-1 

.30 

.4171-3 

.4161-2 

.3441-1 

.3833-1 

.1470-1 

.35 

.3228-3 

.3220-2 

.2593-1 

.2164-1 

.3187-2 

.40 

.2360-3 

.2353-2 

.1839-1 

.1030-1 

-.2373-2 

.45 

.1629-3 

.1624-2 

.1225-1 

.3373-2 

-.4307-2 

.50 

.1057-3 

.1053-2 

.7610-2 

-.3623-3 

-.4345-2 

.55 

.6384-4 

.6357-3 

.4326-2 

-.2015-2 

-.3594-2 

.60 

.3505-4 

.3485-3 

.2160-2 

-.2449-2 

-.2662-2 

.65 

.1657-4 

.1645-3 

.8359-3 

-.2264-2 

-.1834-2 

.70 

.5614-5 

.5548-4 

.1029-3 

-.1833-2 

-.1199-2 

.75 

-.2145-6 

-.2487-5 

-.2464-3 

-.1364-2 

-.7563-3 

.80 

-.2812-5 

-.2826-4 

-.3668-3 

-.9535-3 

-.4643-3 

.85 

-.3547-5 

-.3550-4 

-.3652-3 

-.6342-3 

-.2799-3 

.90 

-.3333-5 

-.3331-4 

-.3092-3 

-.4044-3 

-.1665-3 

.95 

-.2732-5 

-.2729-4 

-.2381-3 

-.2487-3 

-.9794-4 

1.00 

-.2057-5 

-.2052-4 

-.1714-3 

-.1481-3 

-.5689-4 

1.10 

-.9792-6 

-.9765-5 

-.7657-4 

-.4852-4 

-.1855-4 

1.20 

-.3938-6 

-.3925-5 

-.2942-4 

-.1480-4 

-.6189-5 

1.30 

-.1391-6 

-.1386-5 

-.1006-4 

-.4467-5 

1.40 

-.4443-7 

-.4426-6 

-.3169-5 

-.1391-5 

TABLE  8 


PRESSURE  DISTRIBUTION  P  FOR  *  -  -  .1 


T| 

R  -  -  .1 

R  «  -  1 

R  -  -  50 

R  -  -  100 

.00 

-.4733+5 

-.4740+3 

-.5352+1 

-.4555 

-.1867 

.05 

.4220+5 

.4240+3 

.6198+1 

.2102+1 

.2087+1 

.10 

.9162+5 

.9194+3 

.1228+2 

.2974+1 

.2525+1 

.15 

.1119+6 

.1123+4 

.1454+2 

.3004+1 

.2318+1 

.20 

.1126+6 

.1130+4 

.1436+2 

.2684+1 

.1962+1 

.25 

.1017+6 

.1019+4 

.1282+2 

.2272+1 

.1632+1 

.30 

.8495+5 

.8518+3 

.1067+2 

.1885+1 

.1377+1 

.35 

.6671+5 

.6689+3 

.8434+1 

.1572+1 

.1204+1 

.40 

.4958+5 

.4973+3 

.6400+1 

.1342+1 

.1096+1 

.45 

.3497+5 

.3510+3 

.4714+1 

.1187+1 

.1037+1 

.50 

.2340+5 

.2351+3 

.3415+1 

.1089+1 

.1008+1 

.55 

.1479+5 

.1489+3 

.2475+1 

.1033+1 

.9962 

.60 

.8761+4 

.8856+2 

.1835+1 

.1004+1 

.9933 

.65 

.4773+4 

.4867+2 

.1426+1 

.9921 

.9939 

.70 

.2299+4 

.2393+2 

.1182+1 

.9886 

.9956 

.75 

.8739+3 

.9686+1 

.1049+1 

.9893 

.9972 

.80 

.1301+3 

.2259+1 

.9849 

.9915 

.9984 

.85 

-.2016+3 

-.1048+1 

.9608 

.9939 

.9991 

.90 

-.3047+3 

-.2071+1 

.9575 

.9959 

.9996 

.95 

-.2950+3 

-.1969+1 

.9633 

.9975 

.9998 

1.00 

-.2406+3 

-.1421+1 

.9719 

.9985 

.9999 

1.10 

-.1209+3 

-.2206 

.9869 

.9996 

.1000+1 

1.20 

-.4616+2 

.5282 

.9952 

.9999 

.1000+1 

1.30 

-.1193+2 

.8707 

.9987 

.9999 

1.40 

-.1000+1 

.1000+1 

.1000+1 

.1000+1 

TABLE  9 


TANGENTIAL  VELOCITY  V  FOR  k  «  -  .22 


11 

R  -  -  .1 

R  -  -  1 

R  -  -  10 

R  -  -  30 

R  -  -  50 

.00 

.7800 

.7800 

.7800 

.7800 

.7800 

.05 

.6675 

.6675 

.6654 

.6548 

.6670 

.10 

.5573 

.5573 

.5536 

.5356 

.5664 

.15 

.4514 

.4513 

.4467 

.4254 

.4816 

.20 

.3516 

.3516 

.3466 

.3251 

.4106 

.25 

.2595 

.2594 

.2545 

.2347 

.3494 

.30 

.1761 

.1761 

.1715 

.1542 

.2942 

.35 

.1022 

.1022 

.9795-1 

,8329-1 

.2421 

.40 

.3790-1 

.3786-1 

.3414-1 

.2188-1 

.1915 

.45 

-.1690-1 

-.1694-1 

-.2017-1 

-.3038-1 

.1419 

.50 

-.6270-1 

-.6273-1 

-.6547-1 

-.7397-1 

.9361-1 

.55 

-.1002 

-.1002 

-.1025 

-.1095 

.4735-1 

.60 

-.1303 

-.1303 

-.1322 

-.1380 

.3954-2 

.65 

-.1540 

-.1540 

-.1555 

-.1602 

-.3577-1 

.70 

-.1723 

-.1723 

-.1735 

-.1772 

-.7121-1 

.75 

-.1861 

-.1861 

-.1870 

-.1899 

-.1021 

.80 

-.1963 

-.1964 

-.1970 

-.1992 

-.1282 

.85 

-.2038 

-.2038 

-.2043 

-.2059 

-.1499 

on 
•  «  - 

-.2091 

-.2091 

-.2094 

-.2106 

-.1674 

.95 

-.2128 

-.2128 

-.2130 

-.2138 

-.1813 

1.00 

-.2153 

-.2153 

-.2155 

-.2160 

-.1920 

1.10 

-.2181 

-.2181 

-.2182 

-.2184 

-.2061 

1,20 

-.2193 

-.2193 

-.2193 

-.2194 

-.2136 

1.30 

-.2198 

-.2198 

-.2198 

-.2198 

-.2173 

1.40 

-.2199 

-.2199 

-.2199 

-.2199 

-.2190 

TABLE  10 


AXIAL  VELOCITY  G  -  w JI/2CH  FOR  k  -  -  .22 


T| 

R  -  -  .1 

R  -  -  1 

R  -  -  10 

n 

.00 

0 

0 

0 

0 

0 

.05 

.4068-5 

.4067-4 

.3970-3 

.1046-2 

.1767-2 

.10 

.1256-4 

.1255-3 

.1219-2 

.3106-2 

.5238-2 

.15 

.2128-4 

.2127-3 

.2051-2 

.5017-2 

.8373-2 

.20 

.2772-4 

.2771-3 

.2650-2 

.6136-2 

.9946-2 

.25 

.3076-4 

.3074-3 

.2909-2 

.6224-2 

.9368-2 

.30 

.3027-4 

.3024-3 

.2819-2 

.5319-2 

.6522-2 

.35 

.2678-4 

.2676-3 

.2437-2 

.3623-2 

.1613-2 

.40 

.2117-4 

.2114-3 

.1851-2 

.1409-2 

-.4958-2 

.45 

.1435-4 

.1432-3 

.1152-2 

-.1048-2 

-.1269-1 

.50 

.7156-5 

.7125-4 

.4223-3 

-.3514-2 

-.2105-1 

.55 

.2428-6 

.2112-5 

-.2743-3 

-.5813-2 

-.2955-1 

.60 

-.5964-5 

-.5995-4 

-.8966-3 

-.7837-2 

-.3777-1 

.65 

-.1124-4 

-.1127-3 

-.1423-2 

-.9532-2 

-.4540-1 

<>70 

-.1551-4 

-.1554-3 

-.1849-2 

-.1089-1 

-.5222-1 

.75 

-.1883-4 

-.1886-3 

-.2179-2 

-.1195-1 

-.5810-1 

.80 

-.2132-4 

-.2135-3 

-.2426-2 

-.1273-1 

-.6302-1 

.85 

-.2312-4 

-.2315-3 

-.2604-2 

-.1330-1 

-.6701-1 

.90 

-.2438-4 

-.2441-3 

-.2729-2 

-.1370-1 

-.7015-1 

.95 

-.2524-4 

-.2527-3 

-.2813-2 

-.1397-1 

-.7256-1 

1.00 

-.2581 -4 

-.2584-3 

-.2869-2 

-.1414-1 

-.7436-1 

1.10 

-.2640-4 

-.2643-3 

-.2927-2 

-.1433-1 

-.7658-1 

1.20 

-.2661-4 

-.2665-3 

-.2949-2 

-.1440-1 

-.7764-1 

1.30 

-.2669-4 

-.2672-3 

-.2957-2 

-.1443-1 

-.7809-1 

1.40 

-.2672-4 

-.2675-3 

-.2959-2 

-.1443-1 

-.7826-1 

TABLE  11 


RADIAL  VELOCITY  G  -  •  U  FOR  h  «  »  .22 


T1 

R  -  -  .1 

R  -  -  1 

R  -  -  10 

R  -  -  30 

R  -  -  50 

.00 

0 

0 

0 

0 

0 

.05 

.1426-3 

.1425-2 

.1387-1 

.3592-1 

.6070-1 

.10 

.1831-3 

.1830-2 

.1762-1 

.4265-1 

.7140-1 

.15 

.1575-3 

.1574-2 

.1488-1 

.3168-1 

.4998-1 

.20 

.9660-4 

.9650-3 

.8737-2 

.1232-1 

.1101-1 

.25 

.2451-4 

.2442-3 

.1581-2 

-.8637-2 

-.3447-1 

.30 

-.4224-4 

-.4232-3 

-.4968-2 

-.2685-1 

-.7864-1 

.35 

-.9411-4 

-.9418-3 

-.9994-2 

-.4004-1 

-.1164 

.40 

-.1274-3 

-.1275-2 

-.1216-1 

-.4759-1 

-.1447 

.45 

-.1426-3 

-.1426-2 

-.1453-1 

-.4992-1 

-.1626 

.50 

-.1428-3 

-.1428-2 

-.1443-1 

-.4812-1 

-.1702 

.55 

-.1323-3 

-.1323-2 

-.1329-1 

-.4349-1 

-.1685 

.60 

-.1152-3 

-.1152-2 

-.1153-1 

-.3728-1 

-.1595 

.65 

-.9544-4 

-.9544-3 

-.9519-2 

-.3054-1 

-.1450 

.70 

-.7565-4 

-.7564-3 

-.7525-2 

-.2403-1 

-.1272 

.75 

-.5765-4 

-.5764-3 

-.5722-2 

-.1822-1 

-.1080 

.80 

-.4238-4 

-.4238-3 

-.4199-2 

-.1336-1 

-.8887-1 

.85 

-.3014-4 

-.3014-3 

-.2981-2 

-.9486-2 

-.7100-1 

.90 

-.2077-4 

-.2077-3 

-.2051-2 

-.6536-2 

-.5512-1 

.95 

-.1390-4 

-.1390-3 

-.1371-2 

-.4375-2 

-.4163-1 

1.00 

-.9043-5 

-.9042-4 

-.8908-3 

-.2850-2 

-.3060-1 

1.10 

-.3535-5 

-.3534-4 

-.3475-3 

-.1118-2 

-.1526-1 

1.20 

-.1253-5 

-.1253-4 

-.1230-3 

-.3979-3 

-.6840-2 

1.30 

-.4090-6 

-.4089-5 

-.4008-4 

-.1303-3 

-.2730-2 

1.40 

-.1266-6 

-.1265-5 

-.1239-4 

-.4039-4 

-.9338-3 

TABLE  12 


PRESSURE  DISTRIBUTION  P  FOR  r  -  -  .22 


T1 

R  -  -  .1 

R  “  -  1 

R  ■  -  10 

R  -  -  30 

R  -  -  50 

.00 

.1610+7 

.1606+5 

.1341+3 

.6718+1 

.6854 

.05 

.3614+7 

.3606+5 

.2931+3 

.1249+2 

.8849 

.10 

.4213+7 

.4208+5 

.3389+3 

.1371+2 

.9272 

.15 

.3920+7 

.3910+5 

.3126+3 

.1218+2 

.8693 

.20 

.3130+7 

.3122+5 

.2475+3 

.9266+1 

.7502 

.25 

.2157+7 

.2151+5 

.1686+3 

.5927+1 

.5981 

.30 

.1214+7 

.1210+5 

.9252+2 

.2802+1 

.4362 

.35 

.4231+6 

.4215+4 

.2917+2 

.2631 

.2855 

.40 

-.1617+6 

-.1619+4 

-.1734+2 

-.1530+1 

.1642 

.45 

-.5374+6 

-.5365+4 

-.4685+2 

-.2576+1 

.8516-1 

.50 

-.7317+6 

-.7303+4 

-.6170+2 

-.2984+1 

.5467-1 

.55 

-.7870+6 

-.7854+4 

-.6536+2 

-.2912+1 

.7194-1 

.60 

-.7483^ 

-.7466+4 

-.6151+2 

-.2529+1 

.1302 

.65 

-.6553+6 

-.6539+4 

-.5340+2 

-.1987+1 

.2188 

.70 

-.5394+6 

-.5382+4 

-.4357+2 

-.1399+1 

.3254 

.75 

- .4220+6 

-.4210+4 

-.3374+2 

-.8423 

.4385 

.80 

-.3161+6 

-.3153+4 

-.2494+2 

-.3600 

.5485 

.85 

-.2276+6 

-.2270+4 

-.1763+2 

.3166-1 

.64  86 

.90 

-.1582+6 

-.1577+4 

-.1191+2 

.3333 

.7349 

.95 

-.1063+6 

-.1059+4 

-.7657+1 

.5553 

.8060 

1.00 

-.6910+5 

-.6884+3 

-.4619+1 

.7125 

.8623 

1.10 

-.2647+5 

-.2631+3 

-.1146+1 

.8906 

.9368 

1.20 

-.8696+4 

-.8576+2 

.2968 

.9642 

.9752 

1.30 

-.2142+4 

-.2038+2 

.8270 

.9912 

.9928 

1.40 

.1000+1 

.1000+1 

.1000+1 

.1000+1 

.1000+1 

APPENDIX  B 


a.  Pure  Stagnation  Flow, 

Distorted  Stagnation  Flow,  c„  Limiting  Cell  Flow, 
Cell  Flow,  e.  Pure  Wake  Flow. 


STRONG 


Fig.  2:  Critical  Reynolds  number  jR^  versus  velocity  ratio 
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Fig.  3a:  Tangential  velocities  V  versus  dimensionless  variable  T] 

at  Reynolds  number  ! R !  -  1  and  different  velocity  ratios  h. 


G 

Fig.  4b:  Axial  velocities  G  -  vJi/2M  versus  dimensionless  variable  T| 
at  Reynolds  number  I R |  “1  and  different  velocity  ratios  h. 


Pig.  5a:  Radial  velocities  G  ■  -  U  versus  dimensionless  variable 
number  | R I  *  1  and  different  velocity  ratios  h. 
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Fig.  6a:  Pressure  functions  P  versus  dimensionless  variable  T)  at 
Reynolds  number  |R|  “1  and  different  velocity  ratios  h. 


Fig.  7:  Tangential  velocities  V  versus  dimensionless  variable  7] 

at  velocity  ratio  k  =  -  .1  and  different  Reynolds  numbers  R. 
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Fig.  9:  Radial  velocities  G  ■  -  U  versus  dimensionless  variable  T|  at 
velocity  ratio  h  =  -  .1  and  different  Reynolds  numbers  R. 
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Fig.  10:  Pressure  functions  P  versus  dimensionless  variable 
1]  at  velocity  ratio  h  ■  -  .1  and  different  Reynolds 
numbers  R. 
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Fig.  11:  Tangential  velocities  V  versus  dimensionless  variable  T) 
at  velocity  ratio  h  -  -  .22  and  different  Reynolds 
numbers  R. 


.  12:  Axial  velocities  G 
velocity  ratio  h  ■ 


w^W /2fl£  versus  dimensionless  variable  T|  at 
.22  and  different  Reynolds  numbers  R. 
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